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Abstract 

Motivated by concerns regarding possible two-body contributions to the recently-measured in- 
clusive T(nS)^ 7 + X (n=l, 2, 3) direct photon spectra, we report on a new study of exclusive 
radiative decays of these narrow T(nS) resonances into two-body final states W/y, with 1Z a narrow 
resonant hadronic state decaying into four or more charged particles. Such two-body processes 
are not explicitly addressed in the extant theoretical frameworks used to calculate the inclusive 
direct photon spectra, and must also be explicitly inserted into Monte Carlo simulations. Using 
data collected from the CLEO III detector at the Cornell Electron Storage Ring, we present upper 
limits of order 10 -4 for such bottomonium two-body decays as a function of the recoil mass M-r. 

PACS numbers: 13.20.Gd,13.20.-v,13.40.Hq 
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Introduction 



CLEO recently extracted a s from a measurement of the direct photon spectra in 
T(1S,2S,3S) — > ggj^. That extraction was based on a comparison of the ggj width 
to the dominant three-gluon width of the narrow bottomonium resonances. Since the direct 
photon is observable above background only for relatively high energies (E 7 > E hea _ m /2), 
some model dependence is inherent in the determination of the total ggj rate. Given a 
prescription relating the parton- level rate to a s , one can then use that rate to determine a s . 
To extrapolate beyond the experimentally accessible direct photon energy region, CLEO 
relies on theoretical parameterizations of the expected photon energy spectrum in the T 
system 0, [1] to obtain the total direct T — > gg^ decay width relative to the dominant 
T — > ggg width. The theoretical calculations are generally done at the parton level, and 
therefore avoid possible resonant contributions to the photon energy spectrum due to two- 
body decays, e.g., T — > ggj — > jTZ, with 71 some resonant hadronic state. Alternatively, one 
would like to understand hadronization in T — > ggj — > jlZ, which is measured experimen- 
tally as a two-body process. For example, CLEO has recently observed signals in several 
low-multiplicity modes[4j, [5|, and has presented preliminary results consistent with the final 
state T(1S)— > 7 + f/, rf — > n+n^ri, r\ — ► 7r + 7r _ 7r°[6]. We point out, however, that the product 
of (branching ratio) x (efficiency) for all known exclusive modes, combined with the observed 
level of background, implies a signal yield below the statistical sensitivity required to be 
observed as a distinct signal in the inclusive spectrum. 

From the experimental standpoint, the presence of possible digluon resonances opposite 
the photon leads to "bumps" in the otherwise smooth predicted theoretical photon spectra. 
The inability of the current calculations to directly address two-body effects, in part, restricts 
the applicability of T decay models to the region z 7 < 0.92, with z 7 defined as the scaled 
photon energy (= E 1 / E^ eam ) . T decay models are also not reliable above this point because 
one of the emitted gluons is of such low energy that a perturbative calculation cannot be 
trusted jl]. Given that primary glueball candidates are of order 1 GeV in mass, we expect 
the endpoint region of the photon energy spectrum (z 1 > 0.92) to be most susceptible to 
such contamination. Older estimates of ot s based on inclusive radiative photon production 
in T decay using the BLM[9] prescription have consistently yielded values smaller than 
those obtained from different techniques 0]. Recently, it has been realized that color octet 
contributions previously ignored in the older BLM calculation result in estimates of a s in 
excellent agreement with estimates made at the Z-resonance Nevertheless, subtraction 
of possibly-enhanced exclusive contributions to the overall rate would, in principle, result in 
a lowered estimate for the ratio of rates T(T — > ggj)/T(T — ► ggg) and a correspondingly 
larger estimate for a s . To set the scale, given that the typical branching fraction for T — > gg'j 
is of order 10~ 2 , a total resonant enhancement at the level of 5 x 10 -4 would result in a ~5% 
reduction in the estimate of a s at the Z-pole. 

By comparison, a large fraction of all J/ip — > gg^j decays have been identified as two- 



body [7j. Two-body contamination (and larger relativistic corrections) also makes the J/ip 
system somewhat less reliable than the T in estimating a s . This systematic consideration in 
the gg^f analysis motivates our search for radiative decays of the T to resonances: Y(nS)— > 
7?£ (n=l, 2, 3). We concern ourselves with high multiplicity (> 4 charged tracks) final 
states, as we employ the same hadronic event selection cuts in this analysis that we did 
in the gg^y analysis [1]. We note that although two-body branching fractions have been 
observed for, e.g., T(1S)— > 7/2(1270) at the level of 10~ 4 , the fraction of / 2 (1270) decays 
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into >4 charged tracks is only m 3% [3]. 

The analysis, in general terms, proceeds as follows. After selecting a high-quality sample 
of e + e~ annihilations into hadrons using the hadronic event selection cuts of the previous 
analysis [H, we construct the inclusive isolated photon spectra in data taken at both on-Y- 
resonance and off-T-resonance energies (the latter samples are used for systematic checks of 
the overall procedure). A two-body radiative decay of the T will produce a monochromatic 
photon in the lab frame; the energy of the radiated photon E 1 is related to the mass of 
the recoil hadron TZ via Mji = 2E^ eam ^l — z 1 . In the case where the intrinsic width of the 
recoil hadron is much smaller than the experimental photon energy resolution, the measured 
radiative photon energy should be a Gaussian centered at the energy E 7 . For a 1 GeV (4.5 
GeV) recoil photon, this implies a recoil resonance with width typically narrower than 20 
MeV (260 MeV). Not knowing a priori the mass of the hadron TZ, we therefore perform a set 
of fits of the T(nS) photon spectrum to a Gaussian signal, centered at a series of E 1 values, 
and with resolutions corresponding to the known CLEO III electromagnetic calorimeter 
resolution atop smooth polynomial backgrounds, over the range 0.2 < z 7 < l.O. 1 

We construct 95% confidence level upper limits from these fits as a function of recoil mass 
M-R, corrected for the efficiency loss due to the fiducial acceptance of the detector and the 
event and photon-selection cuts that define our data sample. In estimating this correction, 
we assume TZ has spin=0, with a corresponding 1 + cos 2 9 angular distribution for the recoil 
gamma; higher spins will generally give flatter angular distributions. Note that the exact 
form of the efficiency correction due to the event and photon selection cuts varies with the 
decay final states considered for TZ. In the high-momentum region typical of the particles 
recoiling against the direct photon, our per-track charged particle detection efficiency is 
generally >90% over the fiducial volume. To be conservative, we derive our ^-dependent 
efficiency correction from the decay mode yielding the worst reconstruction efficiency. 

This final efficiency-corrected limit is converted into an M^-dependent branching ratio 
upper limit B(-yTZ) by dividing the resulting yield by the calculated total number of resonant 
T events. For the off-resonance running, the distributions are divided by the off-resonance 
luminosity for the sake of comparison between continuum samples. An example of simulated 
signal superimposed on background is given in FiguredJfor the hypothetical process T(4S) — > 
7 + TZ, TZ — > 7r + 7r~7r + 7r~. 

Event Selection 

Event selection criteria in this analysis are identical to those imposed in the previous 
analysis The inclusive photon spectra are therefore identical to those taken from our 
previous analysis as well. The background shape is approximately exponential in the region 



It should be noted that "bumps" in the inclusive photon spectrum can be due not only to resonant two- 
body decays but also to continuum threshold effects such as the crossing of the cc threshold. Given our 
photon energy resolution, all processes of the type e + e~ — > D [ *'D result in photons relatively close 
in energy and produce an apparent enhancement in the region z 7 ~ 0.8. In the previous analysis^, 
we also identified an excess of photons in data as z 7 — > 1 (= £ 7 /£beam)- Further examination of these 
events indicated that they were dominated by continuum production: e + e~ — > 77r + 7r~7r + 7r~, although 
the possibility that the 4-pion state resulted from the decay of an intermediate resonance TZ was not 
investigated. 
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FIG. 1: Scaled photon spectrum for T(4S) — > jTZ, 1Z — > 4ir simulations, for various hypothetical 
1Z masses. The lower (dashed) curve is T(4S) data, while the upper (solid) curve is T(4S) data 
with signal Monte Carlo added on top. The magnitude of i3(T(4S) — > ^ylZ, 1Z — > 4tt) has been 
grossly exaggerated for the sake of presentation (~ 5 x 10~ 3 , well above any observed radiative 
decay fraction into > 4 charged tracks). The six elevations correspond to masses 7.5, 6.5, 4.5, 3.5, 
2.5 and, for the right-most peak, the overlap of a 1Z of mass 1.5 GeV and a 7Z of mass 0.5 GeV, 
respectively. 

0.2 <z^< 1.0. 

Fitting the Inclusive Photon Spectrum 

To extract the possible magnitude of a two-body radiative signal, we step along the 
inclusive photon spectrum over the interval 0.2 < < 1.0, fitting it to a Gaussian with 
width equal to the detector resolution at that value of photon energy, plus a background 
parametrized by a smooth Chebyshyev polynomial. We assume that the intrinsic width 
of the resonance 1Z is considerably smaller than the detector resolution. Our step size is 
determined by the energy resolution of the detector cr E ; we use steps of width <7e/2. 

For photons in the central barrel (| cos# 7 | <0.7, with # 7 the polar angle of the photon 
momentum vector relative to the beam axis) region of the Csl electromagnetic calorimeter, 
the energy resolution over the kinematic interval relevant to this analysis is of order 2%. 
For two-body radiative decays from the T(1S), the photon energy resolution and recoil mass 
resolution function of shown in Figure [2j Curves are similar for the other narrow 

T resonances. 

We use a fourth-order Chebyschev polynomial to describe the background. For each fit, 
the Gaussian fit area A(z 7 ) and fit error <ja{z 7 ) is recorded. Note that the fits are highly 
correlated point-to-point, and that the bin width is much finer than the detector resolution. 
At each step, we use a ±10<te fitting window; the background is expected to be relatively 
smooth over such a limited interval. 
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FIG. 2: Photon energy resolution (solid) and recoil mass resolution for T(1S)— >• "fTZ (dashed). 
Extracting Upper Limits 

To convert the A(z^) distribution obtained from fitting the inclusive photon spectrum 
into a 95% confidence-interval upper limit, we add 1.645 • <ta{zj) point-wise to the A(z^) 
distribution, as a function of photon energy. In this process, since we are interested in 
enhancements in the inclusive photon spectrum, all negative areas from the raw fits are set 
equal to zero, and the corresponding upper limit set to 1.645 • (Ta(z 7 ) at these points. The 
resulting contour for the T(1S) fitting is shown in Figure El The analogous contours for the 
other T spectra look similar. 




FIG. 3: A(z^) + 1.645 • (ta{z^) versus z 7 for fits to the T(1S) inclusive photon spectrum where 
negative points have been mapped to 1.645 • cr^z-y). This plot is the upper limit yield, before 
efficiency correction, computed from statistical errors only. 

We convert the limits as a function of photon energy z 1 into a function of a hypothetical 
resonance recoil mass M-r,. For the purposes of this conversion, the mean values for each 
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running period of -Ebeam are used for each data sample; we neglect the MeV-scale variation 
in beam energies for a particular run period. These values are given in Table I, along with 
the integrated luminosities of the resonance and below-resonance T data samples used in 
this analysis. 



Event Type 


^beam (GeV) 


C (pb- 1 ) 


Resonance T(1S) 


4.73 


1076±11 


Resonance T(2S) 


5.01 


1189±12 


Resonance T(3S) 


5.18 


1228±12 


Resonance T(4S) 


5.29 


6456±65 


Below T(1S) 


4.72 


188±2 


Below T(2S) 


5.00 


396±4 


Below T(3S) 


5.16 


158±2 


Below T(4S) 


5.27 


2278±23 



TABLE I: The mean values of average beam energy (-©beam) an d the integrated luminosity (C) for 
each data sample used in this analysis. 

The resulting M^-dependent contour, for the T(1S), is shown in Figure HI 
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FIG. 4: A(M-ji) + 1.645- a a(Mti) versus M-r for fits to the T(1S) inclusive photon spectrum, where 
negative points have been mapped to 1.645 • ^(M^). Statistical errors included only. 



Efficiency Correction 

We consider two efficiency corrections to the upper limit contour: one due to the fiducial 
acceptance of the detector, and the other due to our event and shower selection cuts. 

For photons in the barrel of the detector, we assume that 1Z is spin zero, which corresponds 
to a 1 + cos 2 6 distribution of the photons in the two-body decays we are considering; higher 
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spins will generally give natter angular distributions [4|. Assuming this angular distribution 
amounts to an « 0.6 uniform angular acceptance efficiency correction factor to our limit. 

In addition to this angular acceptance correction, we assess an efficiency correction due 
to the CLEO III detection efficiency. Not knowing a priori what the decay mode of our 
hypothetical resonance 1Z will be, we have generated 5000-event Monte Carlo samples span- 
ning a wide range of final state multiplicities (all with iV charged > 4) and masses Mr.. In the 
interests of producing a conservative upper limit, we used this study to choose the mode 
with the worst average efficiency. In this manner, we efficiency-correct, as a function of Mr 
mass, our upper limit as a function of z 7 (or E 7 ) before mapping the upper limit into M-r. 
A list of TZ decay modes considered in this analysis and their average efficiencies (averaged 
over the photon energy spectrum from 1.0 GeV < E 1 < 4.5 GeV) is given in Table II. 
We find that the worst efficiency among the decay modes considered was obtained from 
1Z — > 2(K + K~)tc° (Figure E]). We therefore use this efficiency function to determine our 
upper limit contours. 

3031106-023 
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FIG. 5: The efficiency for detecting an T — ► 7 + 1Z, 1Z —* 2{K + K~)ir° event as a function of 
observed photon energy as determined from a 25000-event Monte Carlo sample. Each point 
in this efficiency is obtained from a different 1Z mass hypothesis. This photon energy-dependent 
efficiency correction distribution is used to point-wise correct our upper limit, where the efficiency 
between points in this distribution is estimated by linear interpolation. 
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Results 

To convert the efficiency-corrected upper limit contour into an upper limit on the two- 
body radiative branching ratio B^IZ), we simply divide the efficiency-corrected upper limit 
contour by the total calculated number of T(nS) decays[l|], as shown in Table III. For 
completeness, we also include the results for the T(4S), for which the decay width is expected 
to be nearly saturated by T(4S)— ► BB[7]. The total number of T(4S) events was obtained 
by multiplying the total luminosity of our on- resonance T(4S)-running by the well-known 
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Event Type 


Average Efficiency (e) 


TZ -» K+K--K+-K- 


0.53 


± 


0.03 


72. -> K+K-Vr+TT-TT 


0.53 


± 


0.02 


71 -» K + K-ir + iT-iT ir 


0.54 


± 


0.02 


TZ -> K+K-p+p- 


0.56 


± 


0.02 


1Z -> K + K-p + p-ir° 


0.50 


± 


0.05 


TZ -► K + K-p+p-ir°ir 


0.57 


± 


0.02 


TZ — > p + p _ 7r + 7r _ 


0.62 


± 


0.03 


72 — > p + p~7T + 7T _ 7r 


0.54 


± 


0.05 


72. — > p + p _ 7r + 7r _ 7r 7r 


0.63 


± 


0.02 


72. -» K+K-K+K- 


0.50 


± 


0.02 


TZ -> K + K-K + K--k°tt° 


0.49 


± 


0.02 


TZ — ► p + p~p + p~ 


0.67 


± 


0.02 


TZ — ► p + p~p + p~7T° 


0.65 


± 


0.02 


72. — > p + p _ p + p _ 7r 7r° 


0.63 


± 


0.02 


72. — > 7r + 7r _ 7r + 7r _ 


0.59 


± 


0.02 


72 — ► 7T + 7r _ 7r + 7r _ 7r 


0.65 


± 


0.02 


72 — ► 7T + 7r _ 7r + 7r _ 7r 7r 


0.59 


± 


0.01 


72 — > 7r + 7r _ 7r + 7r~47r 


0.57 


± 


0.02 


72 — > 7r + 7r~7r + 7r~67r° 


0.60 


± 


0.02 


72 — ► 7T + 7r _ 7r + 7r _ 87r° 


0.60 


± 


0.02 


72. K + K-K + K~K + K- 


0.68 


± 


0.04 


72 — > p + p~j> + p~p + p~ 


0.53 


± 


0.04 


72 — > 7r + 7r _ 7r + 7r _ 7r + 7r _ 


0.74 


± 


0.03 


72. -> K+K-K+K-tt 


0.48 


± 


0.02 



TABLE II: Average efficiencies for the reconstruction of various decay modes that could be detected 
in this analysis, obtained by fitting the photon energy-dependent reconstruction efficiencies to a 
straight line in the interval 1.0 GeV < E~ < 4.5 GeV (statistical errors only). Quoted efficiencies 
correspond to 5000-event Monte Carlo samples for which photons are generated with a flat angular 
distribution within the barrel fiducial volume of the calorimeter, and therefore are restricted in 
geometry. The lowest-efficiency final state (K + K~ K + K~ ir° in bold above) is used for setting 
upper limits. 



T(4S) on-resonance cross-section (which we calculate from 111]). The resulting on- resonance 
upper limits B(jlZ) are shown in Figure El 

Given the fact that we have not performed a continuum subtraction on the on-resonance 
inclusive photon spectrum from T decays, it is interesting to compare the structure observed 
in Figure E] with structure observed when we apply the fitting procedure to continuum data. 
Figures El M and [TU] show the resonances' limits of Figure El separately, with the corre- 
sponding continuum distributions overlaid for comparison. We observe a partial correlation 
between the continuum and resonance spectra, suggesting that both spectra have large con- 
tributions from initial state radiation (ISR) photons. 

Applying our fitting procedure directly to the continuum we can obtain limits on the 
cross-section for e + e~ — > 7 + TZ, over the barrel angular acceptance region (Figure ITTT) . 
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FIG. 6: The Af^-dependent B{^1Z) upper limit contours obtained for T — > 7 + 7£, 7£ — >> 4 charged 
tracks for the T(1S), T(2S), T(3S) and T(4S). Limits are obtained by dividing upper limits on 
yield by reconstruction efficiency and number of resonant events, and also incorporating systematic 
uncertainties. All limits are of order B(")1Z) ~ 1CP 4 . 
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FIG. 7: Comparison of the Af-^-dependent Gaussian fit area upper limit A(M-ji) + 1.645- a a(M-ji) for 
the T(1S) versus the below T(1S) continuum; we observe some correlation between the resonance 
and the below-resonance structure. Note that the normalization of the below-resonance upper limit 
curve is arbitrary and has been adjusted so as to allow an easy visual comparison. 
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FIG. 8: Comparison of the M^-dependent Gaussian fit area upper limit A(Mfc) + 1.645 • a^M-ji) 
for the T(2S) versus the below T(2S) continuum; scaling of the continuum as before. 
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FIG. 9: Comparison of the M^-dependent Gaussian fit area upper limit A(Mn) + 1.645 • a^M-ji) 
for the T(3S) versus the below T(3S) continuum; scaling of the continuum as before. 



It is important to note here that a) the angular distribution for continuum initial state 
radiation processes is considerably more forward-peaked than the 1 + cos 2 9 distribution 
we have assumed for the resonance; we have therefore applied a correction based on the 
angular distribution appropriate to ISR, and b) the quantum numbers of particles produced 
in association with ISR photons are different than those produced in radiative decays of 
quarkonium resonances. To set the scale of the continuum cross-section sensitivity, the raw 
ISR cross-section for e + e~ — ► Jjijj + 7 is expected to be ~5 pb in the 10 GeV center-of-mass 
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FIG. 10: Comparison of the Af-^-dependent Gaussian fit area upper limit A(M-ji) + 1.645 • ga(M-ji) 
for the T(4S) versus the below T(4S) continuum; scaling of the continuum as before. 



TABLE III: The total number of calculated T(1S), T(2S), T(3S) and T(4S) events in our data 



regime. Taking into account the efficiency of our event selection requirements and the strong 
forward peaking expected for ISR processes, this corresponds to an expected observed cross- 
section into >4 charged tracks ~ 10~ 4 nb. This value is at the edge of our current statistical 
sensitivity. 

Cross-Check 

In order to check that we are able to identify a signal at a given sensitivity level, we 
embedded pure Monte Carlo signal into data, and performed our fitting procedure on the 
resulting distribution to ensure that we recover the correct signal magnitude in our branching 
ratio upper limit. To do this, hypothetical T(4S)^ 7 + 71, 1Z — > 7r + 7r~7r + 7r~ events were 
embedded into the T(4S) inclusive photon spectrum with branching ratios of the order of 
lCT 5 , 10" 4 , 10" 3 , and 10" 2 under 10 different M n hypotheses: M n = 0.6 GeV, 1.5 GeV, 
2.5 GeV, 3.5 GeV, 4.5 GeV, 6.5 GeV, 7.5 GeV, 8.5 GeV, and 9.5 GeV. The resulting upper 
limit contours derived from applying our procedure to these spectra are show in Figure [121 
We reconstruct all signals within our expected sensitivity (around 10 -4 ) that are within our 



T Resonance 



iV total (T(nS)) (xlO 6 ) 



T(1S) 
T(2S) 
T(3S) 
T(4S) 



20.96 ± 0.06 
8.33 ± 0.04 
5.24 ± 0.06 
6.8 ± 0.2 
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FIG. 11: Mfc-dependent cross-section upper limit contours obtained for e + e~ —* 7 +1Z, TZ — >> 4 
charged tracks for the below T(1S), T(2S), T(3S) and T(4S) continua (nb). This plot is obtained 
by dividing the result of our fitting procedure on the continuum by the off-resonance luminosity. 
The angular correction here is based on the expected distribution appropriate for continuum initial 
state radiation. Systematic errors have also been incorporated into these limits. 

accessible recoil mass range. 
Systematic errors 

We identify and account for systematic errors as follows: 

1. We account for possible systematics in our event and shower reconstruction efficiency 
by using the lowest-efficiency final state considered, and by assuming 1Z has spin=0. 
The angular distributions for spin=0, 1 and 2 two-body decays have been calculated, 
and generally yield natter distributions for higher spins [3]. 

2. We have assessed fitting systematic uncertainties by varying the recoil mass bin width 
(from 20% to 50% of the resolution a) and the order of the background polynomial 
used to parametrize the background (from second-order to fifth-order). Observing 
no statistically significant variation between these extremes, we have used as defaults 
cr=5 bins and a fourth-order background, based on the goodness of fit of the pull 
distributions to a unit Gaussian on the continuum. 

3. For continuum measurements, we assess a uniform 1% degradation of the limit due to 
the luminosity uncertainty as calculated in the previous analysis [D]. 

4. For on-resonance measurements, we degrade the limit uniformly by the uncertainty in 
the calculated number of total resonant events given in Table III. 
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FIG. 12: Upper limit contours derived from applying our procedure to fabricated Monte Carlo 
signal spectra. We reconstruct all input signals withing our sensitivity (~ 10 ) that are within 
our accessible recoil mass range. 

Summary 

As shown in Figure El our sensitivity is of order 10~ 4 across the mass range corresponding 
to 0.2< z 7 <1.0, well above the tabulated branching ratios for any known T — > 7+X, 
/i + /i~/i + /i~+neutrals process. We measure upper limits of: 

B(T(lS)-> 7 + n, n ^> 4 charged tracks) < 1.26 x 10~ 3 , 
£(T(2S)^ 7 + n, n ^> 4 charged tracks) < 9.16 x 10~ 4 , 
£(T(3SH l + n,1Z^>A charged tracks) < 9.69 x 1(T 4 

for all kinematically allowed masses Mn, under the assumption that 1Z is a pseudoscalar. 
Constraining 1.5 GeV < M-r < 5.0 GeV we set a more stringent limit of: 

B(T(lS)-> 7 + K, K ^> 4 charged tracks) < 1.78 x 10~ 4 , 
S(T(2S)^ 7 + 11,11 ^> 4 charged tracks) < 1.95 x 10" 4 , 
£(T(3SH 7 + ^> 4 charged tracks) < 2.20 x 10~ 4 . 

Additionally, we report these upper limits as a function of the mass recoiling against the 
photon, as shown in Figure EJ 

We limit the branching ratio for two-body radiative decays to narrow resonances (<20 
MeV in width) to be < 10~ 4 . We conclude that distortion of the inclusive photon spectrum 
in our previous extraction of ot s due to the possible contribution of such events is negligible. 
The possibility of resonances with widths greater than our experimental resolution has yet to 
be completely addressed. Further work on exclusive multiparticle final states (e.g., 727r + 27r~, 
r y2K + 2K~ , ^K^K and 'yK K ± 7T T ) would help elucidate the nature of such radiative decays. 
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